Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



STIS SPECTROSCOPY OF YOUNG STAR CLUSTERS IN "THE ANTENNAE" GALAXIES (NGC 4038/4039)* 
Bradley C. Whitmore, Diane Gilmore, C. Leitherer, S. Michael Fall, and Rupali Chandar 

Space Telescope Science Institute, Baltimore, MD 21218 



William P. Blair 

Johns Hopkins University, Department of Physics & Astronomy, Baltimore, MD 21218 



FRANgOIS SCHWEIZER 
Carnegie Observatories, 813 Santa Barbara Street, Pasadena, CA 91101-1292 

Qing Zhang 

Cadence Design Systems Inc., San Jose, CA 95134 
AND 

Bryan W. Miller 

AURA/Gemini Observatories, La Serena, Chile 
Draft version February 5, 2008 

ABSTRACT 

Long-slit spectra of several dozen young star clusters have been obtained at three positions in the 
Antennae galaxies with the Space Telescope Imaging Spectrograph (STIS) and its 52"x 0V2 slit. Based 
on Ha emission-line measurements, the average cluster-to-cluster velocity dispersion in seven different 
cluster aggregates ("knots") is <10 km s~F The fact that this upper limit is similar to the velocity 
dispersion of gas in the disks of typical spiral galaxies suggests that the triggering mechanism for the 
formation of young massive compact clusters ( "super star clusters" ) is unlikely to be high velocity 
cloud-cloud collisions. On the other hand, models where preexisting giant molecular clouds in the disks 
of spiral galaxies are triggered into cluster formation are compatible with the observed low velocity 
dispersions. These conclusions are consistent with those reached by Zhang et al. (2001) based on 
comparisons between the positions of the clusters and the velocity and density structure of the nearby 
interstellar medium. We find evidence for systematically lower values of the line ratios [N II] /Ha 
and [S II] /Ha in the bright central regions of some of the knots, relative to their outer regions. This 
suggests that the harder ionizing photons are used up in the regions nearest the clusters, and the 
diffuse ionized gas farther out is photoionized by 'leakage' of the leftover low-energy photons. The low 
values of the [S II] /Ha line ratio, typically [S II] /Ha < 0.4, indicates that the emission regions are 
photoionized rather than shock heated. The absence of evidence for shock-heated gas is an additional 
indication that high velocity cloud-cloud collisions are not playing a major role in the formation of 
the young clusters. 

Subject headings: galaxies: star clusters — galaxies: interactions — galaxies: individual (NGC 4038, 
NGC 4039) 



1. INTRODUCTION 

The discovery of young massive star clusters in merg- 
ing galaxies was one of the important early results from 
the Hubble Space Telescope (Holtzman et al. 1992; Whit- 
more et al. 1993). These clusters have all the attributes 
of young globular clusters, hence providing us with the 
opportunity to study the formation of these systems in 
the local universe. Dozens of other galaxies have now 
been found to harbor similar clusters. Merging galaxies 
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appear to have the largest populations of such "super 
star clusters," though similar clusters are also found in 
smaller numbers in starburst, barred, and even normal 
spiral galaxies (see Whitmore 2003 for a review). This 
discovery has acted as a catalyst for the study of star 
clusters and for the study of star formation in general 
since it now appears that most stars form in clusters 
(e.g., Lada & Lada 2003; Fall et al. 2005). 

What triggers the formation of these young, compact, 
massive clusters? One possibility is that high-velocity 
collisions of the gaseous components in interacting galax- 
ies result in the formation of large-scale shock fronts that 
trigger gravitational collapse and the subsequent forma- 
tion of stars and clusters (e.g., Gunn 1980; Schweizer 
1987; Kang et al. 1990; Olson & Kwan 1990a, 1990b; 
Noguchi 1991; Kumai et al. 1993a, 1993b). Specifically, 
Kumai et al. 1993a suggest that chaotic cloud-cloud col- 
lisions with velocities > 50 - 100 km s _1 may result in the 
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formation of young globular clusters due to local shock 
compression in the Magellanic Clouds. A more recent 
model by Bckki ct al. (2004) suggests that velocities in 
the range 10 - 50 km s _1 are optimal for making clus- 
ters. Recent numerical, smooth-particle-hydrodynamic 
simulations of NGC 4676 ("The Mice") by Barnes (2004) 
have been successful in reproducing the extended nature 
of star formation in this interacting galaxy pair with a 
shock-induced star-formation law. However, it is impor- 
tant to note that much of the extended star formation in 
this model originates in relatively low- velocity shocks due 
to convergent flows, rather than in high- velocity shocks 
created by interpenetrating collisions of galactic disks. 
Some observational support for cluster formation being 
triggered by high-velocity collisions comes from Wilson 
et al. (2000), who find that the strongest CO and IR 
source in the Antennae galaxies coincides with a knot in 
the "Overlap Region" , where two giant molecular clouds 
appear to be colliding with differential velocities of 50 - 
100 km s _1 . 

There is, however, direct evidence against high-speed 
cloud collisions being the dominant triggering mecha- 
nism for the formation of star clusters in the Antennae 
galaxies. Zhang ct al. (2001) made detailed comparisons 
between the positions of star clusters in the Antennae 
and the stellar and interstellar contents of these galax- 
ies, based on observations over a wide range of wave- 
lengths. They found little or no correlation between the 
positions of young clusters (very near their birthplaces) 
and the velocity structure of the surrounding interstellar 
medium, as measured by gradients and dispersions in the 
HI and Ha channel maps at various distances from the 
clusters, as would be expected if cloud motions triggered 
the formation of the clusters (see Figures 15, 16, and 17 
of Zhang et al. 2001). In contrast, Zhang et al. found 
a strong correlation between the positions of the young 
clusters and several other properties, including gaseous 
emission from young stars and the density structure of 
the surrounding interstellar medium as measured by the 
HI and CO intensity maps. The implications of these re- 
sults for the formation of the clusters are discussed fur- 
ther by Fall (2004). 

In the present paper we measure the velocities of sev- 
eral dozen young clusters in the Antennae by using the 
high spatial resolution made possible by the Hubble Space 
Telescope (HST), and from these compute the corre- 
sponding velocity dispersion. While previous ground- 
based studies (e.g., Burbidge & Burbidge 1966; Rubin 
et al. 1970; Amram et al. 1992) mapped out the veloci- 
ties of various "knots" in The Antennae (e.g., labeled A 
through T by Rubin ct al.), the Hubble observations by 
Whitmore & Schweizer (1995) revealed that each of these 
"knots" consists of roughly a dozen or more individual 
star clusters. Long-slit spectroscopic observations with 
the Space Telescope Imaging Spectrograph (STIS) can 
provide a direct measurement of the velocity dispersion 
between clusters within such knots. If clusters formed in 
high-speed cloud collisions, we might expect the resulting 
velocities of the clusters, and hence their dispersion, to be 
large. The main caveat here concerns the possibility that 
some of the kinetic energy of cloud motions is dissipated 
and radiated away during the collisions and hence is not 
fully reflected in the velocity dispersion among the clus- 
ters. If we neglect this dissipation for the moment, the 



observed velocity dispersion among the clusters should 
reflect the velocity dispersion of the clouds from which 
they formed. This is the test we perform here, as a com- 
plement to the one reported by Zhang et al. (2001). 

Throughout the present paper we adopt a Hubble Con- 
stant of Hq = 75 km s -1 Mpc -1 , which for a heliocentric 
systemic velocity of czq — 1630 km s -1 (corresponding 
to 1440 km s -1 relative to the Local Group) places NGC 
4038/4039 at a distance of 19.2 Mpc. At this distance, 
the projected scale is 1" = 93 pc, and the corresponding 
distance modulus is 31.41 mag. We note that Saviane et 
al. (2004) determine a distance of only 13. 8± 1.7 Mpc for 
The Antennae, based on the tip of the red giant branch 
observed in the southern tidal-tail population. This dis- 
tance would imply a peculiar velocity of ^400 km s -1 , 
which would be quite high for a field galaxy in a loose 
group such as NGC 4038/4039. We prefer to use the 
original rcdshift distance for this paper, pending a more 
definitive measurement of the photometric distance. 

2. OBSERVATIONS AND REDUCTIONS 

Long-slit spectra of the Antennae galaxies were ob- 
tained with STIS during the period from 2000 April to 
2001 March, as part of Proposal GO-8170. A variety 
of observations were made, including long-slit Ha ob- 
servations (G750M grating), both long-slit and slitlcss 
observations in the UV (G140L grating), and slitless ob- 
servations in the blue (G430M grating). Most results 
in the current paper are based on the Ha observations. 
More detailed results from observations with the G140L 
and G430M gratings will be reported in a future paper. 
Chandar et al. (2005) also use the G140L spectra in their 
study of the ages of clusters in 27 star-forming galax- 
ies. Their spectroscopically determined ages are in good 
agreement with the photometrically determined ages re- 
ported in the present paper. 

Table 1 presents a variety of measured and derived 
parameters for the cluster candidates discussed in this 
paper. These include the cluster ID number (from Whit- 
more & Zhang 2002), a quality index based on how well 
the Ha and I-band images correspond, absolute mag- 
nitude (My), age estimate (from Whitmore & Zhang 
2002), equivalent width of Ha, velocity, line width, rela- 
tive Ha flux, and [Nil] /Ha and [SII]/Ha line strength ra- 
tios. As discussed below, in cases when the emission ap- 
pears to be diffuse (e.g., the quality = 3 candidates), we 
include the nearest cluster in the table in order to obtain 
information about the physical characteristics within the 
knot (e.g., the mean age of nearby clusters). 

A 52" x 0'.'2 slit was used with the G750M grat- 
ing to cover the spectral region of AA6500 - 7050 A, 
hence including the [Nil] AA6548,6583, Ha, and [SII] 
AA6716,6731 lines. The scale was 0.05" pixel -1 , the 
reciprocal spectral dispersion 0.56 A pixel -1 , and the 
spectral resolution for a nominal point source ~0.9 A 
(FWHM). 

For the observations, three slit positions were selected 
to include as many young clusters as possible, as shown 
in Figure 1. We note the curious linear alignments of 
clusters utilized in visits 6 and 7. It is unclear whether 
these alignments are due to chance, the edges of occult- 
ing dust clouds (which are often found to be filamentary 
and roughly linear), or perhaps shock fronts formed at 
the intersection of the two colliding galactic disks. In 
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any case, the fortuitous alignments allowed us to include 
one to two dozen clusters along each of the three slit 
positions. 

Figures 2a-c show the three long-slit spectra. Note 
that the continuum can only be seen for the brightest 
knots. Figure 3 shows an interesting, peculiar steep- 
gradient feature in the spectrum of Knot B (Visit 6). 
This feature is visible in both the Ha and [Nil] lines. 
The fact that the feature appears tilted indicates that it 
stems from a spatially directed, rather than spherically 
symmetric outflow. There is no obvious corresponding 
feature visible in the Ha and broad-band images. 

The three exact positions of the slit were determined 
post facto by comparing the Ha flux profiles measured 
along the slit with linear brightness extractions from our 
Ha image (Whitmore et al. 1999) taken with the Wide- 
Field Planetary Camera 2 (WFPC2). Figures 4a-c com- 
pare the imaging and spectral profiles for each visit. Us- 
ing this method we were able to determine positions of 
the slit with an accuracy of about O'.'l, good enough to 
deduce which specific objects fell on the slit (Table 1). 
Note that the slightly poorer spatial resolution of the 
WFPC2 Ha image relative to the STIS data makes some 
of the objects appear less prominent in Figs. 4a-c. Fig- 
ures 5a-c show the slit positions superposed on both Ha 
and /-band WFPC2 images. 

Note that not all of the peaks in the Ha profiles are 
due to clusters; some are due to shells and other diffuse 
emission between the clusters. In particular, the frac- 
tion of Ha peaks associated with diffuse emission rather 
than clusters was higher in Knot T, presumably because 
the clusters there tend to be slightly older and have had 
more time to expel their gas. We include a "quality" 
index to represent the degree of certainty of our cluster 
assignments in column 4 of Table 1. This is also why 
we refer to "candidate" clusters rather than one-to-one 
matches. In cases when the emission appears to be dif- 
fuse, we include information about the nearest cluster in 
the table in order to obtain information about the phys- 
ical characteristics within the knot (e.g., the mean age 
of the nearby clusters). Hence it should be kept in mind 
that for the quality = 3 objects, the various properties 
(e.g., the equivalent widths) listed in Table 1 are not 
for the same region that we have sampled from our slit 
spectrum. 

Table 2 summarizes the mean properties of the seven 
primary knots covered by our slit positions. The ages 
of the clusters within these knots were derived using 
the technique described in Whitmore & Zhang (2002), 
but using the Bruzual & Chariot (2003) models rather 
than their 2000 models. Note the correlation between the 
quality index and the age of the clusters. As expected, 
the younger clusters show a better correspondcnse be- 
tween their Ha and continuum morphologies (i.e., lower 
quality index), since the gas has not had much time to 
be spatially displaced from the clusters. We also note 
the extreme youth (~2 Myr) of some of these clusters, 
befitting their large equivalent widths. 

The limited pointing accuracy, and difficulties in deter- 
mining optimal positions for the slit, resulted in typical 
deviations of a few tenths of an arcsecond from some 
desired positions. This caused us to miss some of the 
brightest clusters (e.g., we skirted the edge of Knot B in 
Visit 7). However, with so many clusters in each knot, 



missing one cluster generally meant that we picked up 
another cluster somewhere else along the slit. Hence the 
total number of clusters was roughly the same as origi- 
nally planned. 

The original data were recalibrated using CALSTIS 
with the most recent reference files. Hot pixels were re- 
moved with an IDL routine that replaces deviant pixels 
by the median of the adjacent pixels. Line strengths, 
line widths (FWHM of a Gaussian profile), and veloci- 
ties were measured using the SPECFIT program within 
STSDAS, with results given in Table 1. The simplex al- 
gorithm was used to interactively fit Gaussians to all five 
emission lines simultaneously. Measurements were made 
both in and between the clusters whenever the line flux 
warranted it. Extraction regions ranged from 3 pixels 
(0'.'15) for bright regions to 24 pixels (1'.'2) where the line 
flux was weak. Values for the velocity were taken from 
the Ha peak alone since the weaker lines are absent for 
faint sources. Figures 6a-c present the relative Ha flux 
(solid lines), velocity cz, Ha line width, [Nil] /Ha line 
ratio, and [SII]/Ha line ratio for the three visits. 

A potential complication is the fact that the clusters 
are partially resolved within the ft.'2-wide slit. This 
means that part of the measured velocity is due to the 
off-center position of clusters within the slit. Because 
of this, the velocity dispersions reported below represent 
upper limits. In principle it should be possible to correct 
for this effect by convolving the Ha image with the slit 
transmission function. However, a rough empirical check 
indicated that this effect was barely discernible, whence 
no corrections were applied for the present paper. The 
affect is very small, in any case, as demonstrated by the 
small RMS deviations of radial velocities along the slit, 
as we will discuss later in the paper. 

This complication is more serious for our G140L and 
G430M spectra, most of which were taken with a 0'.'5 
rather than a ft'2 slit. We attempted to circumvent this 
problem by measuring wavelengths for the Milky Way 
components of strong interstellar lines (e.g., C II A1335 
and Si II A1260) in order to determine a wavelength ze- 
ropoint for the spectra. Due to the radial velocity of 
the Antennae, ISM lines arising in the Antennae are red- 
shifted relative to the same lines arising in the Milky 
Way; hence, in principle we should be able to see both 
the Galactic and Antennae components. Figures 21 and 
23 in Whitmore et al. (1999) show that this is indeed 
the case for Knots S and K, that are dominated by two 
of the brightest clusters in the Antennae. Unfortunately, 
the S /N and spatial resolution of the current spectra were 
too low for this method to be applied here successfully. 

A comparison of the newly determined velocities with 
published values for Knots B, C, D, E, F, and T shows a 
relatively large range in absolute differences, especially 
against the early low-dispersion work by Burbidge & 
Burbidge (1966; mean difference of +48 km s _1 , with 
a dispersion of 62 km s _1 ). The agreement is much bet- 
ter for velocities measured by Rubin et al. (1970; mean 
difference of +10 km s _1 , dispersion of 34 km s _1 ) and 
Amram et al. (1992; mean difference of +15 km s _1 , 
dispersion of 18 km s _1 ). A comparison with the HI ve- 
locities from Hibbard et al. (2001) is problematic, due to 
the large mismatch in spatial resolution (i.e., ^O'.'l for 
the HST data and ~10" for the HI data), but a rough 
estimate yields a mean difference of -3 km s _1 and a dis- 
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persion of 48 km s . The differences are in the sense 
that the earlier papers found larger velocities. 

3. RESULTS 

3.1. Cluster-to- Cluster Ha Velocity Dispersions 

A primary goal of this project is to determine whether 
high cloud-cloud velocity dispersions are the triggering 
mechanism for the formation of super star clusters, as 
predicted by some models (e.g., Kumai et al. 1993a, 
Bekki et al. 2004). 

Assuming the gas clouds were already present before 
the collision (i.e., GMCs), one might expect that the re- 
sulting dispersion after the encounter would be compa- 
rable to the original velocity difference between the two 
colliding systems. For example, if portions of two collid- 
ing galaxies encountered each other with relative veloc- 
ities of 100 km s _1 , one might expect a cluster forming 
from the collision of a high-mass cloud of velocity 100 
km s _1 with a low- mass cloud of velocity km s _1 to 
have a resulting velocity near 100 km s _1 . On the other 
hand, if instead the high-mass cloud had a velocity of 
km s _1 , the resulting cluster would have a velocity 
near km s _1 . In practice the situation is likely to be 
much more complicated, with dissipation, multiple colli- 
sions, and turbulence all playing important roles. Hence, 
we would expect the observed cluster-to-cluster velocity 
dispersion to be somewhat lower than the original cloud- 
cloud velocity dispersion. 

A cursory glance at Figures 6a-c reveals immediately 
that the velocity fields along the slit positions are rel- 
atively quiescent, with velocity dispersions ~10 km s _1 
rather than the ~100 km s _1 difference in the orbital 
velocities of the two colliding galaxies. We note that 
these low observed cluster-to- cluster velocity dispersions 
are similar to those found for gas, young stars, and 
young clusters in normal spiral disks (Mihalas & Bin- 
ney 1981). While it is conceivable that the inelasticity 
of cloud-cloud collisions discussed in the previous para- 
graph might reduce the observed velocity dispersion of 
a system of clusters formed by high velocity encounters 
to this degree, it seems too coincidental that the result- 
ing velocity dispersion would essentially match the value 
found in normal disks. 

The low cluster-to-cluster velocity dispersions we have 
observed are more compatible with models where the 
GMCs in the disks of the colliding spiral galaxies are 
the seeds that form the young massive clusters (e.g., Jog 
& Solomon 1992; Elmegreen & Efremov 1997), rather 
than with models where large cloud cloud velocity dis- 
persions are required (e.g., Kumai et al. 1993a, Bekki ct 
al. 2004) . In the former case the cluster-to-cluster veloc- 
ity dispersion would be roughly the same as the original 
cloud-cloud velocity dispersion, i.e., similar to the veloc- 
ity dispersion of the gas in a spiral disk. 

These conclusions are consistent with those reached 
by Zhang et al. (2001), as discussed in the introduction, 
who found little or no correlation between HI and Ha 
velocity dispersions and the presence of young clusters 
in The Antennae. This interpretation is also supported 
by the fact that young massive clusters are distributed 
throughout The Antennae rather than being located only 
in the "Overlap Region," where high-speed collisions are 
most likely to occur. Similarly, young massive clusters 
are also found in quiescent environments such as normal 



spiral galaxies (Larsen & Richtlcr 2000), albeit in smaller 
numbers. 

We have measured the cluster-to-cluster velocity dis- 
persion by selecting knots with three or more measured 
clusters (i.e., Knots B, C, D, E, F, T, and T2), as listed 
in Table 2. In several cases there are velocity gradients 
across the knots which are part of the large-scale veloc- 
ity field (e.g., Knot F in Visit 7). We corrected for such 
gradients by fitting a line to the velocities as a function 
of position in the knot and computing the RMS of the 
velocity residuals. The clusters in a typical knot extend 
~5" along the slit, corresponding to ^500 pc (see Table 
2 for individual knots). 

A potential complication for this procedure is that in 
some cases (primarily Knot C) the velocity field is dom- 
inated by a super bubble surrounding the knot, with the 
central part of the knot appearing blue-shifted by ~50 
km s _1 with respect to its surroundings. This effect can 
be seen most clearly in Figure 2 of Amram et al. (1992), 
where the velocity contours around Knots B and C ap- 
pear to be concentric rings. 

Table 2 shows that the cluster-to-cluster velocity dis- 
persions measured from the ionized gas are low in all 
knots. Averaging over all seven knots yields a mean 
value for the cluster-to-cluster velocity dispersion of 9.6 
km s _1 . Excluding Knot C from the average, because of 
the potential bias introduced by the super bubble, yields 
a negligible difference (i.e., an average value of 9.5 km 
s _1 ). As discussed in §2, the measured velocity disper- 
sions represent upper limits, since no correction has been 
made for the position of each cluster within the slit. 

Globular cluster systems in galaxies show much larger 
velocity dispersions than gas in disks. If the young mas- 
sive clusters in The Antennae are to form a cluster sys- 
tem with kinematics similar to the kinematics of globular 
clusters in typical early-type galaxies, they need to de- 
couple from the gaseous kinematics. This presumably 
happens later in the merger, when the process of vio- 
lent relaxation peaks. The relatively ordered rotational 
velocity field observed in The Antennae by Burbidge & 
Burbidge (1966), Rubin et al. (1970; esp. §IV), and Am- 
ram et al. (1992, esp. Fig. 2) argues that the original 
disks are still relatively intact, and most of the violent 
relaxation lies still ahead. 

The low cluster-to-cluster velocity dispersions we have 
observed in the Antennae beg the question of whether 
the multiple clusters in a given "knot" are gravitation- 
ally bound to each other, and hence may merge together 
to form more massive clusters, as advocated by various 
authors (e.g., Kroupa 1998, Elmegreen et al. 2000, Fell- 
hauer et al. 2005, Bekki et al. 2004). Using the virial 
theorem, and typical values of Mfe not = 4 x 10 6 M Q (each 
knot typically including « 25 individual clusters with 
mass greater than 10 Mq) and Rfcnot = 120 pc for our 
five best defined knots (i.e., B, C, D, F, T), we estimate 
that the binding velocity is « 7 km s _1 (ranging from 5 
to 10 km s _1 for the five knots). While this is roughly at 
the limit where we can make reliable measurements based 
on the current observations, we note that our observed 
upper limit for the cluster-to-cluster velocity dispersion 
is only slightly larger than this estimated binding veloc- 
ity, hence the merging cluster scenario may be viable, at 
least in part, and deserves further study. 

However, it is important to note that the resulting 



Star Clusters in 

structure that forms from the coalescence of all of the 
clusters in a knot would have a radius similar to the orig- 
inal radius of the knot (i.e., w 100 pc), and hence would 
look more like a dwarf galaxy than a globular cluster. It 
may still be possible to form structures with radii similar 
to globular clusters (or the slightly larger "faint fuzzy" 
clusters observed by Brodie & Larsen 2002) if only two 
or three clusters within w 10 pc of each other merge. 

It is also interesting to compare the cluster-to-clustcr 
Ha velocity dispersions with the HI velocity dispersions 
measured by Hibbard et al. (2001) and here listed in the 
last column of Table 2. The measured HI velocity dis- 
persions range from 15 to 30 km s~ 1 , typically a factor 
of two higher than the cluster-to-cluster velocity disper- 
sions. This apparent discrepancy is almost certainly due 
to the poor spatial resolution of the HI maps (~10"). For 
example, Knots T and F both show relatively large Ha 
velocity gradients that we remove when computing the 
cluster-to-cluster velocity dispersion, but that remain in- 
cluded in the high measured HI velocity dispersions. If 
the gradient is not removed the cluster-to-cluster Ha ve- 
locity dispersions are 28.9 km s _1 for Knot T and 31.2 
km s _1 for Knot F; very similar to the values for the HI 
velocity dispersions. 

3.2. Ha Line Widths 

In principle, it should be possible to learn about the 
kinematics of the outflows around clusters and knots 
from the line-width information. Unfortunately, the abil- 
ity to decipher this information is compromised by the 
fact that the 0'.'2 width of the slit is comparable to the 
width of the STIS point spread function. Hence, an 
apparent increase in the line width may be due to the 
presence of diffuse line emission covering the slit more 
uniformly than a point source would, rather than being 
indicative of a kinematic signature. 

According to Table 13.40 of the STIS Instrument 
Handbook, the instrumental FWHM of spectral lines is 
^40 km s _1 for point sources and ~105 km s _1 for very 
extended sources. Note that the latter number is very ap- 
proximate and likely to be a slight underestimate, since it 
assumes a simple rectangular profile without taking into 
account the instrumental wings beyond the edges of the 
slit. In addition, the program used for analysis (IRAF 
task SPECFIT) fits a Gaussian profile rather than a rect- 
angular profile. Hence, direct comparisons with predic- 
tions of the line profile for a uniform diffuse source are 
problematic. 

The measured line widths in the high S /N knots range 
from about 70 to 120 km s _1 (Figure 6). Hence, it seems 
likely that the observed spread is due to the presence of 
a range of objects from partially resolved compact star 
clusters to extended emission from gas that has been ex- 
pelled from the clusters and fills the slit nearly uniformly. 
Because of the fact that the upper limit of the measured 
line profiles is comparable to what is predicted for diffuse 
emission, we conclude that in general there is no evidence 
for outflows with velocities greater than *=s 100 km s _1 . 

Two possible exceptions occur in clusters #7 and #14 
(Visit 6), both with large apparent velocity widths (>150 
km s _1 ) and good signal-to-noise ratios of Ha (i.e., rel- 
ative Ha flux >50 in Table 1). Cluster #14 lies in the 
brightest part of Knot B. Figure 3 shows a peculiar steep- 
gradient feature in its spectrum, with a peak-to-peak ve- 
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locity amplitude of 150 km s _1 . This feature is proba- 
bly responsible for the high velocity width measured in 
Cluster #14. The direct Ha image (Fig. 5b) shows — to 
one side of this cluster — the presence of a large shell-like 
structure that may be associated with the spectral fea- 
ture. The other high-velocity-width feature is seen in 
Cluster #7 (Visit 6), midway between Knots C and D 
(Figs. 2b and 5b). Careful inspection of Figure 2b and 
6b suggests that the superposition of diffuse gas from 
Knots C and D, plus a discrete cluster at a slightly lower 
velocity, may in this case be the explanation for the high 
apparent velocity width. 

Knot T represents a case where diffuse emission proba- 
bly plays the dominant role in determining the line width. 
The correspondence between the Ha and /-band images 
is particularly bad in this knot (i.e., "quality" rating of 
3 for seven out of eight objects in Table 1), probably be- 
cause the clusters are older (mean age ~6 Myr). Hence 
the gas has had more time to escape, resulting in diffuse 
emission which fills the slit more uniformly than a point 
source, resulting in higher apparent line widths (~110 
km s _1 ). 

This interpretation is supported by the fact that the 
Ha and [N II] lines across Knot T (Fig. 2a) show a series 
of wiggles, with peak-to-peak velocity amplitudes of ~60 
km s _1 (see Fig. 6a, around position 500). Note from 
Fig. 6a that the line widths are ~110 km s^ 1 , consis- 
tent with diffuse emission filling the slit. These wiggles 
are relatively rare, and are largely responsible for Knot 
T having one of the largest measured cluster-to-cluster 
velocity dispersions (14 km s -1 , see Table 2). Whitmore 
et al. (1999) inferred similarly high values for the gas 
outflow velocities around Knots S and K based on the 
sizes of the bubbles and age estimates. Hence, it appears 
likely that the observed line wiggles are related to the 
outflow velocities produced by individual cluster winds 
around some of the older clusters. In such cases, the 
stellar cluster-to-cluster velocity dispersions are likely to 
be smaller than the gaseous velocity dispersion. Even so, 
the gaseous velocity dispersion for Knot T is only 14 km 
s _1 . 

The interpretation is simpler for knots that have low 
apparent velocity widths, indicative of cases where the 
emission is coming from compact regions around rela- 
tively point-like clusters rather than from diffuse emis- 
sion that fills the slit and artificially broadens the ob- 
served line width. We can infer that the clusters are 
particularly young in these cases and have therefore not 
had time to blow sizable bubbles of diffuse emission that 
fills the slit. Probably the best example of this is Knot 
D (Fig. 6b), with line widths around 80 km s^ 1 and a 
relatively good correspondence between the Ha and / 
images (i.e., no clear evidence for the presence of shell- 
like structures such as those seen in Knots B and C) . We 
also note that this is one of the younger knots, with a 
mean age for the clusters of ~2.5 Myr according to Ta- 
ble 2. In addition, Fig. 5c shows that the spectrograph 
slit was positioned near the center of this knot, whence 
most of the emission is coming from the bright clusters 
rather than from the diffuse emission in the area. Sim- 
ilarly, Knot F has line widths around 80 km s _1 , and a 
mean age of the clusters of ~1.4 Myr. 

We conclude that the apparent velocities derived from 
line widths are primarily useful for providing informa- 
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tion about whether the emission is relatively point-like 
or diffuse, rather than providing information about the 
actual outflow velocities in the gas. We note that this 
provides a second-order method for estimating cluster 
ages, since the knots with low apparent line widths (i.e., 
more point-like emission) tend to contain the youngest 
clusters. 

3.3. Line-Strength Ratios 

Figures 6a-c show how the [N II] /Ha and [S II] /Ha line 
ratios change along each of the three long-slit positions. 
In each case, what is plotted is the sum of two lines 
(A6584 + A6548 for [N II], A6716 + A6731 for [S II]) 
relative to Ha. 

Several different effects are visible in these figures. 
Overall, the mean [N II]/Ha ratio observed is ~0.4, and 
the mean [S II]/Ha ratio is ^0.2. Faint knots, such as 
seen in the right half of Fig. 6a (Visit 2 - Region T2), 
show a large dispersion due to observational uncertain- 
ties (i.e., clusters with relative Ha flux < 50 in Table 1 
have observational uncertainties ~ 0.25), but are consis- 
tent with these mean values. Secondly, there is evidence 
for systematically lower values of the line ratios in a few 
of the knots where the slit crosses emission of especially 
high surface brightness. This results in the U-shaped line 
ratio patterns in Fig. 6b, Knots B, C, and D. This anti- 
correlation has been seen before in studies of extragalac- 
tic H II regions and diffuse ionized gas (DIG) (Hoopes et 
al. 1999; Greenawalt et al. 1998). Assuming that the DIG 
is photoionized by 'leakage' of leftover ionizing photons 
from star-forming regions, the harder ionizing photons 
must be depleted locally causing a higher mean ioniza- 
tion in the central parts of H II regions (and thus lower 
[S II] /Ha and [N II] /Ha ratios) than in the surrounding, 
lower-surface-brightness diffuse gas. 

The [S II] /Ha ratio is of interest because it has long 
been used as a discriminant for shock-heated versus pho- 
toionized gas (see Blair & Long 1997; Gordon et al. 
1998; and references therein). Usually, a criterion with 
[S II] /Ha = 0.4 as the nominal threshold is used. Shock- 
heated nebulae tend to exhibit higher values of this ratio 
(often considerably higher), while photoionized nebulae 
show lower values (often below 0.2). We see no bright, 
discrete nebulae sampled by the long slits here that ex- 
hibit elevated [S II] /Ha ratios that would be indicative of 
shock heating. Hence, no supernova remnants contain- 
ing bright radiative shocks such as seen in the Cygnus 
Loop in our Galaxy (e.g., Fesen et al. 1985) have been 
detected serendipitously in these observations. 

The absence of obvious kinematic or line ratio evidence 
for shock-heated gas in the vicinity of the young clus- 
ters provides an additional indication that high velocity 
cloud-cloud collisions are not playing a major role in the 
formation of the young clusters. We should note, how- 
ever, that our observations do not rule out the presence 
of low-velocity shocks that heat and compress the molec- 
ular clouds, since cooling takes place primarily through 



molecular lines in the IR in this case, rather than leaving 
a signature in the [S II] /Ha ratio. 

4. SUMMARY 

The Space Telescope Imaging Spectrograph (STIS) 
aboard the Hubble Space Telescope has been used to 
obtain long-slit, high-spatial- resolution spectra at three 
positions in the Antennae galaxies. Observations were 
made in three spectral regions: the far UV (1150 - 1720 
A), the UV-Blue (3800 - 4100 A), and the Red includ- 
ing the Ha, [Nil], and [SII] emission lines. This paper 
focuses on the latter emission lines. The primary results 
are as follows. 

1. The average clustcr-to-cluster Ha velocity disper- 
sion in seven cluster aggregates ("knots") is <10 km s _1 
and, therefore, similar to the velocity dispersion of the 
gas and young stars in the disks of spiral galaxies. 

2. This low cluster-to-cluster velocity dispersion does 
not favor models that require high-velocity cloud-cloud 
collisions to trigger cluster formation (e.g., Olson & 
Kwan 1990a, b; Kumai et al. 1993a, Bekki et al. 2004). 
On the other hand, models where giant molecular clouds 
(GMCs) already present in galactic disks act as the seeds 
from which young star clusters form are compatible with 
the observed low velocity dispersions. This also supports 
earlier findings (Zhang et al. 2001) that there is little or 
no correlation between HI and Ha velocity dispersions 
and the presence of young clusters in The Antennae. 

3. The line ratios of [N II] /Ha and [S II] /Ha across 
knots often show U-shaped patterns, with low values 
where the spectrograph slit crosses the brightest clus- 
ters. This suggests that the harder ionizing photons are 
depleted in regions nearest the clusters, and the diffuse 
ionized gas farther out is photoionized by 'leakage' of the 
leftover low-energy photons. 

4. The observed low values of the [S II] /Ha line ra- 
tio (typically <0.4) indicate that the HII regions associ- 
ated with knots and clusters are photoionized rather than 
shock heated. No supernova remnants containing bright 
radiative shocks have been detected serendipitously by 
our observations. The absence of evidence for shock- 
heated gas in the vicinity of the young clusters provides 
an additional indication that high velocity cloud-cloud 
collisions are not playing a major role in the formation 
of the young clusters. 
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Fig. 1. — Slit positions superposed on the /-band HST image of The Antennae for the three visits using the Ha grating. Capital letters 
identify knots named by Rubin et al. (1970). 
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Fig. 2. — (2a) Spectrum obtained with the Ha grating during Visit 2. Wavelength increases toward the right; the spatial dimension is 
up and down. The lines are (left to right): [Nil] A6548, Ha, [Nil] A6583, and [SII] AA6716, 6731. (2b) Spectrum obtained with the Ha 
grating during Visit 6. The line visible in Knots C and D blueward of the [SII] doublet is He I A6678. (2c) Spectrum obtained with the 
Ha grating during Visit 7. 
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Fig. 3. — Blowup of spectrum of Knot B, including an unusual stcep-gradicnt feature. 
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Fig. 4. — (4a) Comparison of Ha profile along the slit during Visit 2 with an extraction from a simulated slit on the WFPC2 Ha image. 
ID numbers refer to objects in Table 1 and Figure 5a. (4b) Comparison of Ha profile along the slit during Visit 6 with an extraction from 
a simulated slit on the WFPC2 Ha image. ID numbers refer to objects in Table 1 and Figure 5b. (4c) Comparison of Ha profile along the 
slit during Visit 7 with an extraction from a simulated slit on the WFPC2 Ha image. ID numbers refer to objects in Table 1 and Figure 
5c. 
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Fig. 5. — (5a) Slit position and candidate clusters for Visit 2 superposed on the Ha image. Note that the numbers refer to the candidate 
clusters defined in Table 1, and Figure 4. (5b) Slit position and candidate clusters for Visit 2 superposed on the / image. (5c) Slit position 
and candidate clusters for Visit 6 superposed on the Ha image (5d) Slit position and candidate clusters for Visit 6 superposed on the / 
image (5c) Slit position and candidate clusters for Visit 7 superposed on the Ha image (5f) Slit position and candidate clusters for Visit 7 
superposed on the I image. 
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Fig. 6. — (6a) Relative Ha flux (solid lines), velocity cz, Ha line width, and [NII]/Ha and [SII]/Ha line ratios for Visit 2. (6b) Relative 
Ha flux (solid lines), velocity cz, Ha line width, and [Nil] /Ha and [SII]/Ha line ratios for Visit 6. (6c) Relative Ha flux (solid lines), velocity 
cz, Ha line width, and [Nil] /Ha and [SII]/Ha line ratios for Visit 7. 



12 



Whitmore et al. 



TABLE 1 

Candidate Cluster Matches 







Candidate 












Line 


Relative 






Knot 


Object # 


Cluster ID 


Quality a 


M v 


log Age 


2 b log EW(Ha) 


Velocity 


Width 


Ha Flux 


NII/Ha 


SII/Ha 








(mag) 


(yrs) 




(km/s) 


(km/s) 


(km/s) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


VISIT 2 


T 


1 


7025 


3 


—8.1 


7.8 




1623 


82.1 


61 


0.43 


0.17 


T 


2 


6745 


3 


-11.0 


6.7 


2.2 


1614 


102.9 


212 


0.50 


0.17 


T 


3 




3 








1646 


129.5 


198 


0.53 


0.22 


T 


4 


6405 


3 


-11.2 


6.7 


1.9 


1649 


127.0 


187 


0.48 


0.24 


T 


5 


6206 


3 


-10.0 


6.8 


1.9 


1666 


118.6 


99 


0.56 


0.28 


T 


6 


5981 


2 


-12.1 


6.8 


2.2 


1674 


118.6 


270 


0.56 


0.32 


T 


7 


5784 


3 


-11.6 


6.9 


2.1 


1681 


97.1 


290 


0.54 


0.21 


T 


8 


5583 


3 


-10.3 


6.8 


2.4 


1707 


81.3 


281 


0.52 


0.20 


T 
T 


9 
10 












1666 


102.9 


262 


0.49 


0.24 


T2 


11 


4118 


1 


-11.6 


6.0 


3.7 


1699 


198.7 


34 


0.64 


0.77 


T2 
T2 


12 


3635 


2 


-10.0 


6.2 


3.4 


1696 


87.2 


126 


0.40 


0.20 


T2 


13 
14 


3357 


1 


-9.6 


6.5 


3.2 


1676 


127.0 


48 


0.48 


0.57 


T2 


15 


3172 


1 


-12.0 


6.0 


3.6 


1679 


120.1 


38 


0.36 


0.15 


T2 


16 


3118 


2 


-12.9 


6.5 


2.8 


1670 


109.1 


32 


0.34 


0.16 


T2 


17 


2993 


1 


-11.6 


6.5 


3.0 


1694 


123.8 


39 


0.43 


0.51 


T2 


18 


2972 


1 


-11.8 


6.0 


3.7 


1686 


119.1 


49 


0.42 


0.29 


VISIT 6 


D 


1 


2578 


1 


-9.0 


6.0 


3.5 


1446 


75.4 


79 


0.43 


0.27 


D 


2 


2511 


1 


-10.8 


6.5 


3.5 


1458 


81.6 


503 


0.27 


0.11 


D 


3 


2428 


2 


-11.6 


6.6 


3.4 


1477 


70. 


1568 


0.18 


0.07 


D 


4 


2410 


2 


-12.7 


6.8 


2.5 


1479 


92.5 


1304 


0.25 


0.09 


D 


5 


2379 


3 


-11.4 


6.0 


3.5 


1464 


89.1 


1458 


0.21 


0.09 


D 


6 


2314 


2 






3.4 


1420 


93.8 


177 


0.27 


0.15 




7 


2202 


1 


-11.8 


6.0 


3.6 


1489 


156.2 


198 


0.52 


0.29 


C 


8 


2096 


1 


-11.0 


6.0 


3.6 


1466 


85.3 


185 


0.35 


0.24 


C 


9 




3 








1445 


109.7 


626 


0.31 


0.18 


C 


10 


2028 


3 


-12.4 


6.5 


3.2 


1436 


92.7 


1102 


0.28 


0.14 


C 


11 


2002 


2 


-13.3 


6.0 


3.5 


1424 


80.9 


3629 


0.26 


0.05 




12 


1496 


1 


-13.0 


6.0 


3.7 












B 


13 




3 








1469 


116.8 


686 


0.32 


0.12 


B 


14 


1233 


2 


-12.9 


6.5 


3.1 


1444 


152.1 


1490 


0.38 


0.15 


B 


15 


1143 


3 


-12.0 


6.7 


2.1 


1506 


117.9 


241 


0.34 


0.21 


VISIT 7 




1 


9838 


3 


-8.6 


6.9 


2.4 


1648 


94.7 


61 


0.47 


0.07 




2 


9806 


3 


-11. 


6.8 


1.6 


1667 


88.6 


81 


0.25 


0.14 




3 


9745 


1 


-9.9 


6.0 


3.5 


1649 


83.9 


110 


0.39 


0.10 


F 


4 


9196 


3 


-11.2 


6.0 


3.7 


1674 


88.8 


29 


0.49 


0.13 


F 


5 


9152 


1 


-8.3 


6.0 


3.4 


1628 


91.0 


118 


0.48 


0.23 


F 


6 


8996 


1 


-11.9 


6.8 


1.8 


1643 


86.7 


58 


0.35 


0.21 


F 


7 




3 








1627 


71.6 


33 


0.48 


0.21 


F 


8 


8437 


2 


-10.2 


6.0 


3.5 


1610 


71.9 


164 


0.40 


0.19 


F 


9 


8306 


3 


-10.7 


6.0 


3.5 


1622 


82.3 


278 


0.41 


0.14 


F 


10 


8109 


3 


-9.8 


6.0 


3.5 


1617 


82.5 


216 


0.45 


0.23 


F 


11 




3 








1580 


99.5 


215 


0.46 


0.24 


F 


12 


7488 


2 


-11. 


6.4 


3.3 


1580 


88.2 


189 


0.46 


0.23 


F 


13 





3 








1579 


95.2 


161 


0.43 


0.24 


F 


14 


7391 


2 


-11.4 


6.0 


3.7 


1578 


95.3 


76 


0.51 


0.22 


E 


15 


5254 


2 


-10.6 


6.6 


3.2 


1591 


75.5 


124 


0.46 


0.26 


E 


16 


5028 


2 


-12.1 


6.2 


3.4 


1591 


110.8 


232 


0.52 


0.30 


E 


17 


5023 


2 


-11.3 


6.9 


2.5 


1582 


87.8 


239 


0.52 


0.22 


E 


18 




3 








1579 


110.6 


11 


0.30 


0.14 




19 


3884 


1 


-8.4 


6.1 


3.4 


1531 


79.1 


32 


0.37 


0.23 




20 


3193 


1 


-10.9 


6.0 


3.6 


1472 


74.8 


1334 


0.35 


0.22 




21 


3050 


1 


-10.8 


6.0 


3.5 


1514 


79.1 


93 


0.43 


0.15 




22 


2666 


1 


-11.3 


6.9 


2.0 


1457 


94.9 


154 


0.36 


0.24 




23 


2612 


2 


-7.4 


6.8 


2.5 


1474 


97.6 


64 


0.29 


0.22 



a l = Good match between Ha object and cluster; 2 = possible match, 3 = poor match (i.e., Ho diffuse or shell like). Note that in 
cases when the emission appears to be diffuse, we include nearby clusters in the table in order to obtain information about the physical 
characteristics within the knot. 

b Ages based on a slightly modified version of the technique used by Whitmore & Zhang (2002) (i.e., using Bruzualfe Chariot 2003 models). 
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TABLE 2 

Velocity Dispersions and Mean Values of Various Parameters for Selected Knots 





Number of 




Ho Velocity 




Ha Line 






HI Velocity 


Knot 


Clusters 3 


Quality 15 


Dispersions 


log Agc c 


Width 


[Nil] /Ha 


[SII] /Ha 


Dispcrsion d 








(km s- 1 ) 


(yr) 


(km s~ 1 ) 






(km s _1 ) 


T 


10 (7") 


2.9 


14.0 


6.8 


107 (6) 


0.51 (.02) 


0.23 (.02) 


15-30 


T2 


7 (8") 


1.3 


9.8 


6.3 


114 (6) 


0.44 (.04) 




~15 


B 


3 (2") 


2.5 


3.8 


6.6 


117 (1) 


0.36 (.02) 


0.18 (.03) 


~15 


C 


4 (4") 


2.0 


10.5 


6.2 


92 (6) 


0.30 (.03) 


0.14 (.06) 


~15 


D 


6 (2") 


1.8 


15.6 


6.4 


84 (4) 


0.27 (.04) 


0.13 (.04) 




E 


4 (2") 


2.0 


3.1 


6.5 


96 (4) 


0.50 (.02) 


0.26 (.02) 


~30 


F 


11 (10") 


2.1 


10.5 


6.2 


87 (3) 


0.44 (.05) 


0.20 (.01) 


~30 


Mean 


6.4 (5") 


2.1 


9.6 


6.43 


100 


0.40 


0.19 





References. — Numbers in parentheses arc standard deviations of the mean. 
a Numbcrs in parenthesis are the spatial extent of the clusters in the Knot. 
b From Table 1. 

c From Whitmore & Zhang (2004). 
d From Hibbard et al. (2001). 



